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Abstract 



Hadronic events from the data collected with the DELPHI detector at LEP 
within the energy range from 89 GeV to 209 GeV are selected, their jet rates 
are determined and compared to predictions of four different event generators. 
One of them is the recently developed APACIC++ generator which performs 
a massive matrix element calculation matched to a parton shower followed by 
string fragmentation. The four-jet rate is used to measure a s in the next-to- 
leading-order approximation yielding 

a s (M|) = 0.1175 ±0.0030. 

The running of a s determined by using four-jet events has been tested. The 
logarithmic energy slope is measured to be 

dec 1 

= 1.14 ±0.36. 



dlog£ c 



Since the analysis is based on four-jet final states it represents an alternative 
approach to previous DELPHI a s measurements using event shape distribu- 
tions. 
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1 Introduction 

Measurements of hadronic multijet rates in electron-positron annihilation provide an 
excellent test of perturbative quantum chromodynamics (QCD). They can be confronted 
with predictions of QCD-based hadronisation models and allow a precise determination 
of the strong coupling a s . Furthermore, the study of multijet production originating 
from QCD processes is essential for the understanding of the background to four-quark 
production in W + W~ or ZZ decays and also for understanding the background in the 
search for new phenomena. Here we report the final measurements of 2-, 3-, 4-, and 5-jet 
rates using all data collected by DELPHI during the years 1993 to 2000. The 4-jet rate 
is used to determine a s . 

Until 1995 the large electron-positron storage ring LEP at CERN operated at centre- 
of-mass energies around the Z resonance. Due to the high cross-section the total number 
of hadronic events collected during this part of the LEP1 phase is about 2.5 million. 
Analysing LEP1 data enables precise measurements of the strong coupling and detailed 
comparisons of different methods for extracting a s , see e.g. [1]. From autumn 1995 
onwards the centre-of-mass energy was continuously increased and finally reached about 
209 GeV in October 2000. During the LEP2 programme DELPHI collected a total of 
about 12000 hadronic qq events at centre-of-mass energies between 130 GeV and 209 
GeV. The statistics of hadronic events collected at LEP2, though small compared to that 
gathered near the Z resonance, are sufficient for the measurement of jet rates and for a 
determination of the strong coupling a s , see e.g. [2]. Analysing both LEP1 and LEP2 
data gives access to the energy dependence, the running of the strong coupling and thus 
to a direct test of asymptotic freedom. 

In Sec. 2 the selection of hadronic events, the reconstruction of the centre-of-mass 
energy, the correction procedures applied to the data and the suppression of W + W~ and 
ZZ events (and other four-fermion background) are briefly discussed. Sec. 3 presents 
the applied jet clustering algorithms, the measured jet rates and the comparison of the 
data with predictions from hadronic event generators. In Sec. 4, the measurement of a s 
based on the 4-jet data is presented. As in all analyses using topological information from 
hadronic events, the error on the value of a s is dominated by theoretical uncertainties. 
Here we determine a s by applying second order perturbation theory with an optimised 
renormalisation scale. In Sec. 4, the a s measurements along with studies of different 
choices of the renormalisation scale and the investigation of the running of a s from LEP1 
and LEP2 data are presented. 

2 Selection and correction of hadronic data 

The analysis uses data taken with the DELPHI detector at centre-of-mass energies 
between 89 GeV and 209 GeV divided into 14 energy bins. Data entering the analysis, 
including the integrated luminosities collected at these energies and the cross-sections of 
the contributing processes, are summarised in Table 1. 

DELPHI is a hermetic detector with a solenoidal magnetic field of 1.2T. The tracking 
detectors in the barrel part (starting from the beam pipe) are a silicon micro-vertex de- 
tector (VD), a combined jet/proportional chamber inner detector (ID), a time projection 
chamber (TPC) as the main tracking device, and a streamer tube detector (OD) in the 
barrel region. The forward region is covered by silicon mini-strip and pixel detectors 
(VFT) and by the drift chamber detectors (FCA and FCB). 
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Table 1: Data entering the analysis: the columns show the mean centre-of-mass energies 
E cm , the years of data taking, the integrated luminosities, the cross-sections 
for qq (before and after the cut on the effective centre-of-mass energy yjs' rcc > 
0.9 ■ E cm , described in the text), W + W~ , and neutral boson pair production 
(from Zhtter6.21 [3]) and the number of selected qq events after the cuts 
described in the text. 
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Table 2: Criteria for track- and event selection, p is the momentum, Ap its error, r 
the radial distance to the beam-axis, z the distance to the beam interaction 
point (LP.) along the beam-axis, <p the azimuthal angle, -/Verged the number 
of charged particles, thrust the polar angle of the thrust axis with respect to 
the beam, E to t the total energy carried by all measured particles, y/s' TCC the 
effective centre-of-mass energy, E cm = 2E beam = ^Js the nominal centre-of- 
mass energy, and D 2 the discrimination variable, defined in Eq. 1. The first 
two cuts apply to charged and neutral particles, while the other track selection 
cuts apply only to charged particles. 

The electromagnetic calorimeters are the high density projection chamber (HPC) in 
the barrel, and the lead-glass calorimeter (FEMC) in the forward region. The hadron 
calorimeter (HCAL) is a sampling gas detector incorporated in the magnet yoke. Detailed 
information about the design and performance of DELPHI can be found in [4,5]. 

In order to select well measured charged particle tracks, the cuts given in the upper 
part of Table 2 have been applied. The cuts in the lower part of the table are used to 
select e + e~ — > Z/^y — > qq events and to suppress background processes such as two- 
photon interactions, beam-gas and beam-wall interactions, leptonic final states, and, for 
the LEP2 analysis, initial state radiation (ISR) and four-fermion (4f) background. 

At energies above 91.2 GeV, the large cross-section of the Z resonance peak raises 
the possibility of hard intial state radiation (ISR) allowing the creation of a nearly on- 
shell Z boson. These "radiative return events" constitute a large fraction of all hadronic 
events. The ISR photons are typically aligned along the beam direction and usually 
escape detection. In order to evaluate the effective hadronic centre-of-mass energy ^fs 1 
of an event, considering ISR, an algorithm called SPRIME is used [6]. SPRIME is based 
on a fit imposing four-momentum conservation to measured jet four-momenta (including 
estimates of their errors). The hypotheses of single and multi photon radiation are tested 
based on the x 2 obtained in the corresponding constrained fits. 

Figure 1 shows the 189 and 200 GeV effective centre-of-mass energy spectra as com- 
puted with SPRIME for simulated and measured events passing all but the yj s' rcc cut. 
A cut on the reconstructed centre-of-mass energy yj s' rcc > 0.9 • E cm is applied to discard 
radiative return events (see Table 2). Two-photon events, leptonic events and events due 
to leptonic or semileptonic W + W~ decays are strongly suppressed by the cuts. The re- 
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189 GeV 




Vs rec [GeV] 
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Figure 1: Reconstructed centre of mass energy ^s' rec before all cuts except the one on 
\/s' rcc compared to QCD and four-fermion simulations. The small differently 
shaded areas in the bottom right of the plots indicate the size of WW and 
neutral boson pair background. 



maining background from these types of events was found to be negligible in the following 
analysis. 

Since the topological signatures of QCD four-jet events and hadronic 4f events are 
similar, no highly efficient separation of the two classes of events is possible. Furthermore 
any 4f rejection implies a bias to the shape distributions of QCD events, which needs to 
be corrected with simulation. In this analysis a cut in the discrimination variable D 2 [7] 
is applied to separate four-jet QCD events from hadronic W + W~ decays. All events 
are forced into a four-jet configuration by a clustering algorithm. From the resulting 
four-momenta of the pseudo-particles the following quantity is calculated: 

D 2 = min {(Mij - M w ) 2 + (M kl - M w ) 2 } (1) 
= (12; 34), (13; 24), (14; 23). 

The discrimination variable D 2 is based on a comparison of invariant dijet masses to the 
nominal mass of the W boson. The minimum difference for all possible jet pairings (ij, kl) 
is expected to be small for events arising from boson pair production. Figure 2 shows 
the distribution of D 2 at 205 and 207 GeV, compared to the simulation of contributing 
processes. Events from W + W~ or neutral boson pair production cluster at small values 
of D 2 , while e + e~ — > Z/7 — > qq events extend to larger values of D 2 . Demanding D 2 > 
900 GeV 2 leads to an efficient suppression of W + W~ and neutral boson pair events. All 
remaining 4f contributions are estimated by using Monte Carlo generators and subtracted 
from the measurement. The simulations are normalised using the cross-sections given 
in Table 1. The quoted (Jww values correspond to a W mass of 80.35 GeV. For the 
simulation of WW and ZZ events the following generators were used: 

• EXCALIBUR [8] generates four-fermion final states through all possible electroweak 
four-fermion processes. The generator includes the width of the W and Z bosons. 
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Figure 2: Distribution of D 2 for accepted data at 207 GeV before applying four-fermion 
cuts, compared to qq and 4f simulation. 



QED initial state corrections are implemented using a structure function formal- 
ism [9]. EXCALIBUR also includes a Coulomb correction [11] for the CC03 WW 
production [10]. 

• PYTHIA 5.7 [12] is a general-purpose Monte Carlo generator for multi-particle pro- 
duction in high energy physics. As a general-purpose generator it does not contain 
the detailed modelling of all the specific corrections that are contained in the dedi- 
cated four-fermion generators. 

For the central result the EXCALIBUR generator was applied while the difference be- 
tween PYTHIA and EXCALIBUR was used to estimate the systematic uncertainty on 
the background subtraction (see Sec. 4.5). 

Detector and cut effects are unfolded with simulation. The influence of detector effects 
was studied by passing generated events (JETSET/PYTHIA [13] using the DELPHI 
tuning described in [14]) through a full detector simulation (DELSIM [4]). These Monte 
Carlo events are processed with the reconstruction program applying selection cuts as 
for the real data. In order to correct for cuts, detector and ISR effects, a bin-by-bin 
acceptance correction C, obtained from e + e~ — > Z/j — > qq simulation, is applied to the 
data: 

Kfi) gen,noISR / \ 

Wei.i — TJ~T\ ' \ Z > 

"\Ji) acc 

where /i(/i) g en,noisR represents bin i of the shape distribution / generated with the tuned 
generator. The subscript noISR indicates that only events without large ISR (-y/i — 
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V s rcc < 0.1 GeV) enter the distribution. /i(/i) acc represents the accepted distribution / 
as obtained with the full detector simulation. 



3 Jet rates 

Jet clustering algorithms are applied to cluster the large number of particles of a 
hadronic event into a small number of jets, reflecting the structure of hard partons of the 
event. Most clustering algorithms in e + e~ annihilation apply a recursive scheme based 
on an ordering variable d^, a distance measure and a merging scheme indicated by © 
in the following, all being functions of the four-momenta p of two objects i and j. The 
algorithms start with a table of particles representing the initial objects. The pair of 
objects with the smallest dy is considered for merging. These two objects are merged 
into one new object by applying the merging scheme Pk = Pi ® Pj, provided that the 
distance measure is smaller than some given maximum separation y cut . This step is 
repeated with the two particles % and j replaced by the combined object k. After each 
iteration the ordering variables d^ have to be recalculated. The iteration stops if only 
one object remains or if all distance measures are larger than y cut . 

The remaining objects are called jets and the number of jets n is a function of the 
cutoff parameter y cut . 

The n-jet rate, R n (y C ut) gives the fraction of n-jet events relative to all events. By 
definition: 

«,(»»,.) = 1. (3) 



The details of the clustering algorithms used in this analysis are defined below. 

3.1 JADE 

The JADE algorithm [15] is based on the same distance measure and ordering variable: 

_ 2E i E r (l-cos9 ij ) 
u>ij — yij — 2 ' I v 

E vis being the visible energy, which would be the centre of mass energy E cm for a perfect 
detector, E iy Ej being the energy of the objects i and j and 9^ being the angle between 
Pi and pj. 

The merging scheme simply adds the four momenta of Pi and pf 

Pk=Pi® Pj =Pi+ Pj ■ (5) 

There are shortcomings within JADE arising from the choice of the distance measure 
y^. For events with soft gluons radiated off the quark and the antiquark, there are 
kinematical regions where JADE combines the soft gluons first. The resulting "phantom" 
jet has a resultant momentum at large angle to the initial quarks and may point to a 
region where no particles exist. 

3.2 DURHAM 

In case of the DURHAM or k± algorithm [16] the distance measure dij and the ordering 
variable y^ are the same but they are now changed from mass to normalised transverse 
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momentum k i . 



2-min{E 2 ,Ef}-(l-cos%) 
dij = y„ = ^ • ( 6 ) 



vis 



This choice mitigates the shortcomings of the JADE algorithm. 



3.3 CAMBRIDGE 

The CAMBRIDGE algorithm [17] is a modified k± clustering algorithm similar to 
the DURHAM algorithm. It is designed to preserve the advantages of DURHAM while 
reducing non-perturbative corrections at small y and providing better resolution of jet 
substructure. CAMBRIDGE is based on the same distance measure yij as DURHAM 
(Eq.5, 6). The ordering variable dij is a function of the angle between the objects i and 
j ( "angular ordering" ) : 

c^ = 2-(l-cos%) . (7) 

If y%j > 2/cut the object with the smaller energy is stored as a jet and deleted from the 
event table ("soft freezing"). If y^ < y cut the objects are merged into a new object. The 
iteration stops if only one object remains or if all distance measures y^ are larger than 

Ucut- 



3.4 Results 

Figure 3 shows the CAMBRIDGE four-jet rate as a function of y cut from the 
207 GeV data before and after the D 2 > 900 GeV 2 cut and underlines the separation 
power of D 2 . The data are found to be in good agreement with the simulation. 

The remaining amount of four-fermion background is subtracted to obtain the final 
data points given in Figures 4, 5 and 6 showing the jet rates i?2, -R3, R& and R$ as 
determined with the JADE, DURHAM and CAMBRIDGE jet algorithms at 91 GeV and 
for a sample of LEP2 energies. Within errors, the 2-, 3-, and 4-jet rates show a good overall 
agreement at all energies with the generator predictions tuned to data at the Z resonance. 
Figures 7 and 8 show a detailed comparison between CAMBRIDGE jet rates and Monte 
Carlo predictions. Several models, tuned to DELPHI data, are available [14, 18] and are 
used within this analysis: 

• PYTHIA 6.1 is a parton-shower model with explicit angular ordering, followed by 
string fragmentation [12]. 

• ARIADNE 4.08 performs a colour-dipole cascade, followed by string fragmenta- 
tion [19]. 

• HERWIG 6.1 is a coherent parton-shower model, followed by cluster fragmenta- 
tion [20]. 

• APACIC++ performs a massive leading-order (LO) matrix element (ME) calcula- 
tion for 3-, 4-, and 5-jet final states, matched to a parton-shower and followed by 
Lund string fragmentation [21-23]. The APACIC++ parameters have been tuned 
to DELPHI data measured at the Z resonance. 

The precise LEP1 data in particular allow a critical judgment of the precision of tuned 
Monte Carlo models [18]. The parton-shower model PYTHIA tends to overestimate the 
3-jet rate and to underestimate the 4-jet rate at large y cut , see also Figure 4. To cure 
the lack of multijet events a calculation of the underlying matrix elements has been 
performed. APACIC++ also tends to overestimate the 3-jet rate but predicts more 4-jet 
events at small y cut . By taking quark mass effects into account the agreement with the 
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data improves somewhat. The parton-shower generator HERWIG also gives an acceptable 
agreement with the data. The best overall agreement is obtained with the colour-dipole 
model ARIADNE. At LEP2 energies the deviations are obscured by the larger statistical 
errors. Within errors all models show good agreement with the data. Note that the 
errors shown are statistical only and that neighbouring bins are correlated. Considering 
the experimental errors and model uncertainties, no significant excess of multijet events 
at higher energies is observed. 



207 GeV 




Figure 3: Four-jet rate (R4), determined with the CAMBRIDGE algorithm, from raw 
data at 207 GeV as a function of y cu t, compared to the simulation of the 
contributing processes. Top: before cuts against four-fermion background. 
Bottom: After cutting at D 2 > 900 GeV 2 . 
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91 GeV 




Figure 4: Jet rates (R ) at 91 GeV as a function of y cut compared to the prediction of 
PYTHIA 6.1 



12 



91 GeV 



U 1.2 
1 



DELPHI 
APACIC++ v.1.05 M =0 
APACIC++ v.1.05 M X) 
PYTHIA 6.1 
ARIADNE 4.08 
HER WIG 6. 1 




a) 




O 1.2 
1 



91 GeV 





• DELPHI 



APACIC++ v.1.05 M q =0 
APACIC++ v.1.05 M q >0 



PYTHIA 6.1 
ARIADNE 4.08 
HERWIG 6.1 



b) 



gj 1.4 

a 

1 

0.8 




Figure 7: Jet rates determined with the CAMBRIDGE algorithm, a) 3- jet rate at 
91 GeV. b ) 4-jet rate at 91 GeV. The upper inset shows the corrections C& et 
applied to the data. The central plot shows the jet rates with their statistical 
error in comparison with different Monte Carlo predictions. The lower inset 
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Figure 8: 4-jet rate at 189 GeV determined with the CAMBRIDGE algorithm. The 
upper inset shows the corrections det applied to the data. The central plot 
shows the jet rates with their statistical error in comparison with different 
Monte Carlo predictions. The grey area in the central plot shows the already 
subtracted background of WW and ZZ events. The lower inset shows the 
jet rates normalised to the data. The band indicates the statistical and 
systematical uncertainty of the data. 
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4 Determination of ct s 

The strong coupling constant, a s , is determined from the four-jet rate, by fitting an 
a s -dependent QCD prediction folded with a hadronisation correction to the data. 



4.1 NLO predictions 

The next-to-leading-order (NLO) expression of the four-jet rate is given by: 

R A (y) = B(y) ■ a] + [C(y) + 2B(y) ■ b ■ ln^] • a 3 s + . . . . (8) 

Where x M = fi 2 /Q 2 , fi being the renormalisation scale, Q the centre-of-mass energy of the 
event, b = (33 — 2ny)/127r and rif the number of active flavours. Eq. 8 shows the explicit 
dependence of R 4 on the renormalisation scale /i. The coefficients B(y) and C(y) for the 
DURHAM and the CAMBRIDGE algorithms are obtained by integrating the massless 
matrix elements for e + e~ annihilations into four-parton final states, performed by the 
NLO generator DEBRECEN [24,25]. The i? 4 results obtained with the JADE algorithm 
are not used for the a s determination because of phantom jets and larger hadronisation 
corrections. 

Figure 9 shows the dependence of R4 on x^. For small values of the overshoot of 
the NLO expression changes into an underestimate of the observed/measured R4. Thus 
small values of x^ are expected when fitting the data, suggesting important contributions 
of higher-order corrections of 0(a^). 




Figure 9: Predictions (left for DURHAM, right for CAMBRIDGE) of the four-jet rate 
i?4 at 91 GeV using DEBRECEN and Eq. 8 for various values of x^ at fixed 
a s = 0.118. Illustrated is the change of the prediction with varying x^. The 
analytical calculations are compared to the parton level prediction using the 
PYTHIA generator. 
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Figure 10: Distribution of the hadronisation corrections to the four-jet rate. The plots 
show the ratio of the four-jet rates after and before simulation of the frag- 
mentation, evaluated with different Monte Carlo models. 



4.2 Hadronisation 

Before comparing Eq. 8 with the data and fitting its parameters a s (and the 
transition of coloured partons into colourless hadrons has to be accounted for. This 
transition has been simulated using Monte Carlo fragmentation models. For each centre- 
of-mass energy the QCD prediction is multiplied by the hadronisation correction 

r (TP \ — ^had'(-^cm) / q\ 

^hatH-C/cmJ — rSim./p \ ' W 
J part [.-^cm) 

where f^'(E cm ) (fp^l'(E cm )) is the model prediction on the hadron (parton) level at the 
centre-of-mass energy E cm . The parton level is defined as the final state of the parton 
shower created by the Monte Carlo event generation. 

The matching of ME calculations with a parton shower within APACIC++ allows the 
tuning, performed at LEP1 energies, to be extrapolated. Thus APACIC++ is the only 
ME generator available at LEP2 energies. Therefore APACIC++ is taken as the reference 
model. The scatter of results in a s , when using different Monte Carlo generators is added 
to the systematic error. Figure 10 shows the ratios between hadron level and parton level 
as a function of y cut from different generators. Using the CAMBRIDGE algorithm the 

ratio Rhadrcm J gparton & we&ker 

y cut dependence than the same ratio determined by 

using DURHAM. 



4.3 Dependence on the renormalisation scale ii 

The explicit dependence of a s derived from Eq. 8 on the renormalisation scale arises 
from the truncation of the perturbative series after a fixed number of orders. Within 
perturbative QCD is an arbitrary parameter. A conventional scale setting called 
"physical scale" is the choice x fl = \. However, several other proposals for evaluating the 
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renormalisation scale are available in the literature. Two of them are investigated within 
this analysis: 

• Method of effective charges (ECH) [26]: 

In 0(ot z s ) perturbation theory, the ECH scale value has to be chosen in such a way 
that the third-order term vanishes: 

C(y) + 2B(y)-b -\nx fl = 0. (fO) 

• Principle of minimal sensitivity (PMS) [27]: 

The PMS optimisation amounts to the determination of the renormalisation scale 
value, which minimises the sensitivity of the theoretical prediction with respect to 
its variation: 

-?-[C(y) + 2B(y)-b -\nx„]=0. (11) 

Within both theoretical scale-setting methods the scale function of y cut . The 

uncertainty of the scale is conventionally estimated by a scale variation within an ad hoc 
chosen range. 

In perturbative QCD the x^ dependence of the prediction for an observable R would 
vanish in the all orders limit only. It has been shown in [1] that an excellent description 
of precise mz data can be obtained by fitting simultaneously a s and x^. In the same way 
a simultaneous fit of a s and x^ to the jet rates was performed to account for the missing 
higher-order calculations. The fitted scale is called the experimentally optimised scale 
rr° pt . The results of the scale-setting methods are shown in Figure 11. Experimentally 
optimised scales for different fit ranges (indicated by the error bars in y cut -direction) and 
for several hadronisation models are compared with the theoretical scale evaluations. The 
fit ranges for x° pt are varied between y cut = 0.05 and y cut = 0.0005. Below y cut ~ 0.0005 
the perturbative expansion is expected to become invalid, above y cut — 0.05 the number 
of events entering becomes too small to perform the fit. For y cut > 0.001 small values 
of Xn are preferred and for y cut near 0.01 the theoretical scale evaluations are of the same 
magnitude as the experimentally optimised scales. 
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Figure 11: Optimised renormalisation scales: The lines show the y cu t dependence for 
theoretically optimised scales at E cm = 91 GeV. The dots give results for 
experimentally optimised scales. The error bars in the horizontal direction 
indicate the fit range, the different symbols represent hadronisation correc- 
tions applied by different Monte Carlo models. 
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4.4 Fits to LEP1 data and a precise measurement of q. s (M^) 

As discussed above, for each measurement of a s the renormalisation scale has to be 
chosen. To determine the experimentally optimised scale a two-parameter fit of Eq. 8 
with a s and x° pt as free parameters is performed to the four-jet rate. Table 3 shows the 
results for x° pt . 

fx 



algorithm 


fit range 


T opt 
M 


DURHAM 


0.001 - 0.01 


0.015 


CAMBRIDGE 


0.001-0.01 


0.042 



Table 3: Experimentally optimised scales. 



Figure 12 shows results of fits using Eq. 8 and the fit ranges given in Table 3 for 
DURHAM and CAMBRIDGE and for both physical (x M = 1) and experimentally opti- 
mised scales. While the fit with experimentally optimised scales results in a good agree- 
ment with the data over two orders of magnitude in y cut , the fit results with physical 
scale show a y cut dependence inconsistent with the measurement. 

Figure 13 presents results of the a s fits as a function of y cnt for different scale evaluation 
methods. For fits with physical scale the resulting a s values show a strong dependence 
on the choice of y cut . Within the investigated range a s varies from about 0.1 to 0.13. 
Theoretically optimised scales (ECH and PMS) improve the situation, but for small values 
of ?/cut where the theoretical scales increase, a s shows again a strong dependence on y cnt . 
Choosing experimentally optimised scales cures the problem. With this choice the a s 
results are independent of y cut , and furthermore results for DURHAM and CAMBRIDGE 
are in good agreement. Experimentally optimised scales are therefore considered as an 
accurate tool to perform a consistent measurement of the strong coupling from four-jet 
rates. 

The jet rate data, as shown, for instance, in Figure 4, are highly correlated. Therefore 
a second fit is performed to just one single bin in y cut with a s as the only free parameter 
using the fixed scales of Table 3. This final fit is performed at y cnt = 0.0063 for both the 
DURHAM and CAMBRIDGE algorithms. As shown in Figure 13 the fit results are very 
stable in the vicinity of this y cut value. 

The total error on a s (M|) is estimated by considering the following experimental and 
theoretical uncertainties: 

• Variations of the track and event cuts given in Table 2: N chaxged > (7 ± 1), E tot > 
(0.50 ± 0.05) • E cm and (25° ± 5°) < thrust < (155° ± 5°). 

• In order to account for a remaining dependence on y cut , the working point is varied 
in the range 0.0016 < y cut < 0.01. 

• The difference between fit results in a s when exchanging the hadronisation model is 
considered as an estimate of the error due to simulation: This error already includes 
quark mass effects since the APACIC+- 1- model takes b-quark masses into account. 

• To estimate the theoretical error the scale is varied around its optimised value: 
0.5 ■ x° pt < < 2 • x° pt as in [1], covering the scatter of experimentally optimised 
scales obtained with different fit ranges and for different hadronisation models, see 
Figure 11. 

• While b-quark mass effects are included in the hadronisation corrections performed 
with APACIC++ the DEBRECEN generator used to compute the coefficient func- 
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tion in Equation 8 is available only for the massless case. From recent investigations 
of b quark mass effects [18] it has been evaluated that these can shift the result by 
as much as 1.8%. Conservatively a contribution to the uncertainty of this size has 
been added. 

The statistical error, the uncertainties obtained by varying track and event cuts and 
by varying y cut are combined into the experimental error. Table 4 summarises the con- 
tributions to the error on the a s (M z ) measurement and Table 5 contains the a s (M z ) 
results. Within the experimental error the results obtained by using the DURHAM or 
the CAMBRIDGE algorithm are consistent. The total error on the measurement is 3.0% 
for DURHAM and 2.6% for CAMBRIDGE. If the scale is varied around its optimised 
value within the larger range 0.25 • x° pt < < 4 • x° pt the contribution to the error on 
a s (M z ) due to the x^ variation has to be increased from 0.0014 to 0.0085 for DURHAM 
and from 0.0007 to 0.0037 for CAMBRIDGE. 



contribution to error 


DURHAM 


CAMBRIDGE 


statistical error 


0.00045 


0.00050 


cut variations 


0.00041 


0.00020 


working point variation 


0.0011 


0.0008 


total experimental error 


0.0012 


0.0010 


MC model exchange 


0.0023 


0.0017 


b mass effect 


0.0021 


0.0021 


total had. error 


0.0031 


0.0027 


x M variation 


0.0014 


0.0007 


total error on a s (M z ) 


0.0036 


0.0030 



Table 4: Contribution to the error on a s (M z ) for DURHAM and CAMBRIDGE. 



observable 




± 


exp. 


± 


hadr. 


± 


scale 


DURHAM 


0.1178 


± 


0.0012 


± 


0.0031 


± 


0.0014 


CAMBRIDGE 


0.1175 


± 


0.0010 


± 


0.0027 


± 


0.0007 



Table 5: Results in a s {M 2 z ) for DURHAM and CAMBRIDGE 



20 




Figure 12: Fits to the four- jet rate R 4 measured at the Z resonance using different 
scale evaluation methods. Top: the distributions. The hatched curve shows 
the results for the experimentally optimised scales. Bottom: the ratio 
i?4 ata /i?4 t . The grey bands show fit results with the physical scale (x^ = 1), 
the cross-hatched bands for experimentally optimised scales (x^). 
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Figure 13: Dependence of a s on y cut : The grey bands give fit results with physical scale 
(x^ = I), the lines with theoretically optimised scales (ECH, PMS) and the 
rectangles with experimentally optimised scales (x ^). 
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4.5 Measurement of the running of a s 

To investigate the energy dependence of the strong coupling constant a s , the fits to the 
four-jet rates (with optimised scales, as determined in Sec. 4.3) are repeated at all centre- 
of-mass energies listed in Table I. In order to account for the lower statistics of the LEP2 
data, the working points are shifted to smaller values of y cut which, however, are still in 
the range of stable a s results: DURHAM y cnt = 0.0040, CAMBRIDGE y cnt = 0.0025. 

The determination of a s at different energies allows the predicted scale dependence of 
the coupling due to higher order effects to be tested. Starting from the renormalisation 
group equation: 

E 2 cm ^~ = PM = -a 2 s (b + ha s + ...), (12) 
cli^ cm 

the logarithmic energy slope is obtained: 

rU da ;' S = ~- 2 P(as) = 26 (l + ^a s + - ..) = 2b (l+ ^ + . . . V (13) 

dlog(£ cm ) a* b V % log (EgJA 2 ) J 

with b = (33 — 2nf)/l2n, bi = (153 — 19n/)/247r 2 . In leading order the logarithmic 
derivative (Eq. 13) is independent of a s and E cm and twice the coefficient b of the f3 
function (2&o = 1.22 for rif = 5). Evaluating the equation in second order results in a 
small dependence of the derivative on a s . Thus A QCD and an appropriate energy scale 
have to be chosen in order to calculate a single value of the logarithmic derivative which 
can be compared with the experimental result. Using E cm = 150 ± 60 GeV (the average 
energy of our measurements), A = 230 MeV(corresponding to a s (M z ) = 0.118), and 
71/ = 5 one obtains da~ l /d \ogE cm = 1.27 ±0.10. 

The experimental value of daJ 1 /dlog£' cm as obtained from fitting the function 
l/(a\ogE cm + b) to the measured a s values is in good agreement with the QCD ex- 
pectation (Table 6 and Figure 14). 

The following contributions to the systematic error on the logarithmic derivative are 
considered: 

• Since the acceptance corrections C acc (Eq. 2) are correlated between all energies, a 
possible systematic error would have only a reduced influence on the logarithmic 
energy slope of a s . Still the acceptance correction is energy- dependent. To evaluate 
the corresponding systematic error, the difference between the correction factor C 
at the Z pole and at LEP2 energies is added to C at the three energies near the Z 
resonance at 89, 91, and 93 GeV and the fit is repeated. The full difference in the 
slopes found with or without this change is considered as the contribution to the 
systematic error of the logarithmic slope due to the acceptance correction. 

• At LEP2 energies the cut in the reconstructed centre-of-mass energy is changed from 
\/ s rec — 0-9 ' E CTa to \fs' vcc > 0.8 • E cm and the fit is repeated. The difference in the 
logarithmic slopes is taken as the contribution to the systematic error. 

• The treatment of 4f background is an important source of systematic uncertainties. 

— The 4f simulation is performed by using alternatively PYTHIA or EXCALIBUR, 
the full difference being included as the systematic error. 

— For the subtraction of 4f background the cross-section is varied by its total error 
on ±1.5%. 

— The cut in the discriminating variable is varied: D 2 > (900 ± 100) GeV 2 . 

• The renormalisation scale is varied at all energies: 1/2 • x° pt < < 2 • a:° pt . 
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Effects due to track and event selections are regarded as fully correlated between 
the energies and thus neglected. Table 6 contains the results of the da^ 1 /dlogE CTa 
measurements for the DURHAM and CAMBRIDGE algorithms and the corresponding 
statistical and systematical errors. 



Observable 


da^/dhgEcm 


DURHAM 
CAMBRIDGE 


1.21 ±0.26 ±0.20 
1.14 ±0.25 ±0.26 


QCD expectation 


1.27 



Table 6: Results of the da~ 1 /dlogE cm measurements for the DURHAM and CAM- 
BRIDGE algorithms. The theoretical expectation is calculated in second or- 
der. 



5 Summary 

Hadronic jet rates in electron-positron annihilation have been measured by DELPHI at 
centre-of-mass energies between 89.4 and 209 GeV. The data agree with the expectation 
from QCD-based event generators. No indication of a significant excess of multijet events 
at high energies is found. 

The strong coupling constant has been determined from the four-jet rate in 0(a 3 s ). A 
variety of methods to solve the renormalisation scale problem has been investigated. A 
consistent measurement of a s can be performed by using experimentally optimised scales. 
The results obtained with two different jet-clustering algorithms agree. The final result 
quoted is obtained by applying the CAMBRIDGE algorithm, since this algorithm has 
small third-order contributions, and shows a smaller dependence on the renormalisation 
scale: 

a s {M 2 z ) = 0.1175 ±0.0030 (tot). (14) 

The result in 0(a 3 ) is statistically uncorrelated and in good agreement with previous 
DELPHI measurements [1] and also with the world average value [28]. The a s result is 
also in good agreement with recent a s measurements of the OPAL [29] and ALEPH [30] 
collaborations based on four-jet rates measured at the Z resonance using 0(o? s ) calcu- 
lations combined with the resummation of large logarithms. The scale-setting methods 
obtained in [1] are confirmed. 

The comparison of a s as measured at the Z and at higher energies confirms that the 
energy dependence of the strong coupling is consistent with QCD expectation. Results 
from DURHAM and CAMBRIDGE are consistent. The logarithmic energy slope, again 
obtained from CAMBRIDGE and again statistically uncorrelated to the result of the 
0(a 2 s ) analysis presented in [31], is measured to be 

-^£- = 1.14 ±0.36 (tot), (15) 



while the QCD prediction for this quantity is 1.27. The measurement is in good agreement 
with previous measurements in 0(a 2 ) [2,31]. 
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Figure 14: Energy dependence of a s as obtained from jet rates with experimentally 
optimised scales. The errors shown are statistical only. The band shows 
the QCD expectation when extrapolating the world average [28] to other 
energies. The dashed lines show the result of the 1 / log y/s fits. 
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